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Abstract 

In complex systems such as turbulent flows and financial markets, the dynamics in long and short 
time-lags, signaled by Gaussian and fat-tailed statistics, respectively, calls for a unified description. 
To address this issue we analyze a real dataset, namely, price fluctuations, in a wide range of tem- 
poral scales to embrace both regimes. By means of Kramers-Moyal (KM) coefficients evaluated 
from empirical time series, we obtain the evolution equation for the probability density function 
(PDF) of price returns. We also present consistent asymptotic solutions for the timescale depen- 
dent equation that emerges from the empirical analysis. From these solutions, new relationships 
connecting PDF characteristics, such as tail exponents, to parameters of KM coefficients arise. The 
results reveal a dynamical path that leads from Gaussian to fat-tailed statistics, furnishing insights 
on other complex systems where akin crossover is observed. 

PACS numbers: 05.10.Gg, 05.40.-a, 89.65.Gh 
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One of the main problems in statistical physics consists in the study of macroscopic 
changes of systems in which fluctuations play a central role, e.g., diffusion and noise-induced 
transitions. The Fokker-Planck equation (FPE) [l| provides a powerful tool for dealing with 
such problems and has been used in many different fields in natural sciences, including solid- 
state and plasma physics, quantum optics, chemical and nuclear reaction kinetics, molecular 
biology and population dynamics. Financial data have also been described as stochastic 
processes governed by Langevin or FPEs. These efforts are of utter relevance due to the 
strongly complex fluctuating dynamics of financial time series, which poses new challenges 
to model the dynamical laws responsible for the observed statistical properties. In that 
direction, for example, the anomalous diffusive properties and second-order correlations of 
orice fluctuations (e.g., see Q]), have been addressed through multivariate [3] and non-linear 
4| models. However, these attempts are usually built phenomenologically and, as far as 
they are manifold, entail drawbacks for unmasking the processes that rule the underlying 
dynamics. 

A typical feature of complex systems is the existence of non-trivial structures on different 
timescales. In particular, the dynamics of price fluctuations in long and short timescales, 
signaled by Gaussian and fat-tailed probability density functions (PDFs), respectively, has 
been usually treated in the literature separately, and a unified description in the full range 
of timescales is still lacking. 

Our present goal is two-fold. The first one is to grasp, through a non-parametric method, 
the underlying stochastic dynamics of price fluctuations in order to unveil the driving mech- 
anisms within a unified framework. The second one is to determine analytically, from the 
stochastic equations that emerge from the first step, the family of PDFs that encompasses 
the observed timescale dependent ones. To these ends, we obtain the evolution equation 
for the PDF of price returns through the estimation of Kramers-Moyal (KM) expansion 
coefficients. We follow the work by Friedrich and collaborators ji], [f| , which exploits a cor- 
respondence between financial market dynamics and hydrodynamic turbulence [71] assuming 
the existence of a flux of information towards finer scales. This approach has been applied 

i n n n 

before for developed markets, although for limited time windows [a, [8|, [9|, [10|. Here, as a 
relevant example, we scrutinize the daily and intraday time series of Ibovespa, the financial 
index of the Brazilian stock market, which is not fully understood, despite typifying ma- 
jor emergent markets. We address a large hierarchy of time-lags, ranging from months to 
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minutes. From the measured KM coefficients, we are able to reproduce the full evolution of 
empirical histograms of price returns, embracing the crossover from Gaussian to strongly fat- 
tailed PDFs when going from large to short timescales. We also present consistent solutions 
of the resulting FPE. They belong to the class of generalized Student-t distributions (also 
known in recent literature as g-Gaussians [111]) that comprises the non- stationary invariant 
PDFs observed in both asymptotic limits. 

We investigate the timescale evolution of the PDF of logarithmic price increments (re- 
turns) Ax. We consider, as a general evolution equation for those PDFs, the KM expansion 
of a master equation, valid for Markov processes [lj]: 



with D^ k > = M^ k \Ax, r, At)/ Ar/k\, being M^ k ' the moments of the conditional PDFs, i.e., 
= / dAx'(Ax' - Ax) k P(Ax',r + At\Ax,t). Following the insight provided by cascade 
models in turbulence, we consider a logarithmic time scale defined as r = log 2 (Ato/Ai), 
where At is a reference time-lag. 

For the analysis of Ibovespa, we select three datasets: 3960 deflated daily closing prices, 
in the term 02Jan.1991-28Dec.2006, 37984 15-minute cotes in 21Jan.1998-31Mar.2003 and 
794310 30-second cotes in 01Nov.2002-19Jul.2006. The timescale At (r = 0) is set as 32 
(trading) days. In what follows, the measured returns are given in units of the standard 
deviation 032 of the respective data sample at 32-day time-lag. 



Markovianity was investigated by evaluating the Chapman-Kolmogorov equation [1]. We 
found that it holds, thus validating our approach. The KM coefficients D^ k \Ax, r) were esti- 
mated directly from data series by means of their statistical definition, given by Eq. (J2j) . Con- 
ditional PDFs P(Ax%, T2IAX1, n) = P(Ax2,T2] Axi,T\)/ P(Ax\,ti), with r 2 > T\, were ob- 
tained from the data sets by building the histograms for the joint PDFs P(Ax 2 , t 2 ; Ax\, Ti), 
computed over pairs of returns Ax,, incident at the same initial time. The first coefficients 
]j( k ) W ere computed for r = (r x + r 2 )/2 and At = t 2 — t\. For each couple of values (r, At), 
we found that L>« and function of Ax, follow, in very good approximation, linear 




(1) 



where the coefficients D^(Ax, t) are defined as 



D {k \Ax,T) = hm D w {Ax,t,At) 



(2) 
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FIG. 1: Coefficients D<$ (a) and (b) ciS ct function of Ax, evaluated at At 2 = 109 min and 
Aii = 125 min (hence r = 7.9 and At = 0.2). The abscissa axes are rescaled by the standard 
deviation 032 of the returns at 32-day time-lag. Solid lines correspond to the intervals used for fits. 

and quadratic laws, respectively, as illustrated in Fig. [TJ Namely, 

Z) (1) = -hAx + h , (3) 
D {2) = b 2 [Ax} 2 + b 1 Ax + b . 

Similar behaviors for D^ l > and have been observed for linear and logarithmic increments 



of indexes and exchange rates involving US, Germany and Japan markets 0. 



3,Q- By fitting 

the linear and quadratic ansatz to the data, we obtained the parameters {ai,bj} for each 
(r, At). For fixed r, the limit At — > in Eq. ([2]) was achieved by extrapolation of the 
parameters {di,bj} as a function of Ar. For the fourth order KM coefficients as 
a function of Ax, fourth order polynomial fits were performed, the limit Ar — > being 
consistent with vanishingly small D^ 4 ^(Ax,r). Therefore, according to Pawula theorem lj 
the KM expansion ([1]) can be truncated after the second order, thus reducing to the form of a 
FPE. The limiting values {oj, bj} determine the Ax-dependence of drift, D^ l \ and diffusion, 
D^ 2 \ coefficients. The r-dependence of the limiting parameters is presented in Figs. [2J131 

Parameters a® and b\ describe the deviation of the PDFs from symmetry around zero. 
Indeed, this feature is observed for larger scales, where both parameters present non-null 
values, although with large fluctuations. In contrast, for intraday scales, and b\ are 
negligible. Drift parameter a\ remains approximately constant along daily and intraday 
time-scales. Diffusion parameter 6 2 increases with r from near zero up to a limiting value, 



signaling a crossover in the underlying dynamics at the monthly scale. Meanwhile, diffusion 
parameter bo presents a sustained exponential decay as r increases. Similar behavior of 60 
was also reported for exchange rates ^] , although for a shorter range of timescales. Despite 
the exponential-like decay of bo, it can not be neglected at large r: representing the additive 
noise component, it provides stability to the stochastic process for small Ax. In fact, the 
variance a 2 , also follows an exponential-like decay 2 _Tr , as shown in Fig. [3](c), thus setting 
the reference level for 60 ■ Two limiting regimes associated with slightly different slopes are 
observed for both quantities (see Fig.[3]^c)), suggesting that b (r) is related to a 2 .. In the limit 
of small t both quantities are characterized by 7 = 1, corresponding to normal diffusion, 
while in the limit of large r, 7 >1, showing the onset of superdiffusion in the high-frequency 
regime. For r > 11, data restrictions for obtaining small At results prevent the estimation 
of the parameters. 

The resulting FPE explicitly reads: 



d T P = d Ax ( ai P) + d 2 Ax ((b + hAx + b 2 (Ax) 2 )P) , 



(4) 



where the r-dependence of the parameters was smoothed according to the ansatz plotted in 



Figs. 



Eq. (jl]) was numerically integrated by means of a FTCS scheme 
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FIG. 2: Dependence of drift coefficient parameters a\ and clq on timescale r. Analysis was per- 
formed with daily (circles), 15-min (squares) and 30-sec (triangles) data series. 



5 




FIG. 3: Dependence of diffusion coefficient parameters 62,^1 and 60 on timescale r. Symbols as 
in Fig. [2j In (c), for comparison, we include <t 2 (t) = a^/a^ 2 (small symbols) together with its 
asymptotic behaviors (thin lines) for large and short timescales as predicted by Eq. ([7]), with 7 =1 
and 1.17, respectively. 

fit to the empirical histogram at At=128 days (r = —2), was used as initial condition. The 
evolution was carried out down to the scale of 30 seconds, the highest time resolution of 
our data. For r >11, further evolution of the FPE, was performed by extrapolation of the 
r-dependence of the coefficients. In Fig. HJ we show the PDFs of returns (rescaled by a T ) 
generated by the FPE, together with the empirical ones. Their agreement is remarkably 
good, in the full range of data, strongly supporting our estimation of KM coefficients. 

Within Langevin dynamics, D^> and are related to the deterministic and random 
forces, respectively [l]. Despite the wide range of analyzed time-lags, the intensity of the 
harmonic restoring force, given by D^(Ax,r) ~ — ai(r)Ax, remains almost constant for 
t > —2. This means that the relaxation mechanisms that are governed, among other factors, 
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by constraints, flux of information, stock liquidity, and risk aversion, act similarly at diverse 
temporal scales. 

The evolution of the diffusion coefficient presents more distinctive features. For most 
timescales, D^(Ax,t) ~ &o( r ) + ^(rM^x] 2 is dominated by the state-independent and 
quadratic components, associated to additive and multiplicative noises, respectively. Due to 
the cumulative character of the fluctuations, the additive component 60 increases with At. 
Meanwhile, the change of 62 to a higher plateau for small At indicates large multiplicative 
effects in that region, that fade away in the opposite limit of large time-lags, as expected. 
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FIG. 4: PDFs of normalized returns. Comparison between the numerical solution of FPE (JH) (solid 
lines) and the empirical data (symbols). Timescales correspond to r = —2 to r ~ 15, from bottom 
to top. Price returns are scaled by a T and PDFs shifted, for better visualization. A Gaussian fitted 
to the data for r = —2 (At=128 days) is the initial condition. 



7 



This means that the endogenous behavior of the market, which spontaneously creates the 
amplification response mechanism to price fluctuations, presents different typical levels for 
micro and macro timescales. 

The presence of multiplicative noise is known to be an ubiquitous mechanism to gen- 
erate fat-tailed PDFs as steady-state solutions {13I . . It was found out that for a large 
set of control parameters, power-law tails prevail, with exponent depending on the ratio 
(Z1/&2 while being bo independent. Furthermore, Langevin dynamics mapped onto stochastic 
multiplicative processes with reinjection gives rise to power-law PDFs P(Ax,r) ~ 

from which the expression /i = 01/(262] was derived [if]]. From these previous results, 
the observed plateaux of a\ and 62 suggest a r-invariant tail of the PDF. In fact, the FPE 
(J1J) admits invariant solutions as discussed below. 

Assuming an exponential law for bo(r), steady values of a\ and 62, and neglecting bi, the 
solution of FPE (0} is 

P(Ax,t) ~ l/(6 (r) + 6 2 [Aa;] 2 )^ +1)/2 , with (5) 
ix = 1 + - B/2)/b 2 , (6) 

where B = —UJbp does not depend on time. Notice that Eq. differs from the expression 
reported in [16]. The solution ((SJ) is a g-Gaussian (with q = 1 + 2/(/i + 1)) [ll| . In the 
particular case b 2 — > 0, the solutions become of the Gaussian form. 
The rescaled variance <J 2 {r) = cr 2 /cr 2 2 * s 

a 2 (r) =b (r)/(a 1 -B/2-b 2 ). (7) 

In the limit of large time-lags, the diffusion coefficient is dominated by the state inde- 
pendent term b that obeys b (r) = 6 (0)2~ r . Substitution of the numerical values of b (0) 
and ai, into Eq. (JTj) yields, in very good approximation, a 2 (r) = 2~ r (normal diffusion in 
the linear time scale) in agreement with numerical results (see Fig. [2(c)). As a consequence, 
in that limit, the evolution equation recovers Gaussianity, ruled by a balance between the 
deterministic harmonic force and the time dependent additive noise. In the opposite limit 
of large r, b 2 attains a non-null steady value, while bo(r) ~ 2 -Tr with 7 > 1. Also in this 
case, Eq. (JTJ) predicts an asymptotic behavior of a 2 in agreement with empirical values, as 
depicted in Fig. 13(c). 
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FIG. 5: Dependence of exponent \x on timescale r. Symbols correspond to fitting of Eq. ([5]) to the 
empirical histograms, the dotted line to the asymptotic value given by Eq. ([6]). 

As the evolution of the parameters other than 6 is slow, especially for r > 1, we also 
checked if the PDFs could be effectively described by the ansatz (j3J) in an extended temporal 
regime. The outcomes of the fits match well the respective empirical PDFs, for almost any 
r. For t ~ 0, the significantly non-null value of bi, imposes a correction to the g-Gaussian 
form, yielding asymmetry. The values of \i resulting from (least squares) g-Gaussian fits and 
are displayed in Fig. [5] together with the asymptotic value obtained from Eq. (]6]). For small 
r, the increasing value of the effective fi points the onset of the Gaussian regime. For large 
r, the empirical tail exponents tend to a steady value in good accord with the theoretical 
ones. 

Power- 



regimes 



|4, 



aw tails are often quoted in the literature for financial assets in high-frequency 
, • Ansatz of the g-Gaussian form have been proposed before for high-frequency 



18j | and for daily [9] logarithmic returns from a phenomenological approach. Meanwhile, 
in our case, they arise naturally from the evolution equation obtained empirically from the 
evaluation of KM coefficients throughout time-lags. 

Let us recall that the results for r > 11 (At >4 min), in Fig. HI were obtained by 
extrapolation of the r-dependence of the coefficients. However, a good foresight of the 
empirical histograms arises from FPE evolution up to r=14.75 (At=l min). On one hand, 
this ensures the reliability of the estimated parameters; on the other, it accounts for the 
predictability of intraday statistics, due to the well known existence of memory effects in the 
high-frequency regime of price returns. It is worth to mention that the observed deviation 



of the empirical histogram for the smal 



beginning of a non-Mar kovian regime 19] . 



est analyzed time-lag (30 seconds) expresses the 
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In sum, we have disclosed the evolutionary pattern of the empirical PDFs of price returns, 
from Gaussian to long-tailed regimes. Beyond the system under study, the scope of our 
results may be appropriate for a larger financial scenario. Moreover our findings are not 
restricted to the complexity of financial data, but m ay p rovide insights for systems on 



general physical contexts, such as in turbulent flows 20], where similar Gaussian to q- 
Gaussian crossovers have been observed. 
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